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Abstract 

The interplay between disorder and compressibility in Ising mag- 
nets is studied. Contrary to pure systems in which a weak compress- 
ibility drives the transition first order, we find from a renormalization 
group analysis that it has no effect on disordered systems which keep 
undergoing continuous transition with rigid random-bond Ising model 
critical exponents. The mean field calculation exhibits a dilution- 
dependent tricritical point beyond which, at stronger compressibility 
the transition is first order. The different behavior of XY and Heisen- 
berg magnets is discussed. 



*Laboratoire associe au C.N.R.S. (U.A. n° 17) et aux Universite Paris VII et Paris VI. 
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1 Introduction 



The effects of tlie magneto-elastic coupling (MEC), between the spin and the 
lattice vibrational degrees of freedom, have been investigated extensively [1- 
10]. Basic thermodynamics implies that for magnetic systems with a diverg- 
ing specific heat the presence of the MEC will cause an instability. Therefore 
the divergence is averted by a preemptive first-order transition. Mean-field 
theory (MFT) and solutions of simplified models [3-6] have vindicated this 
conclusion. In particular MFT predicts the existence of a tricritical point 
separating second- and first- order transitions, as the intensity of the MEC is 
increased [5] . Other types of behavior were predicted upon changing bound- 
ary and other conditions. However, they will not be discussed here [11]. An 
even better understanding was achieved with the application of the renormal- 
ization group (RG) approach to these systems [7-9]. Calculations in d—A-e 
dimensions, have shown the behavior to depend crucially on the number of 
components m of the order-parameter. For example the ferromagnetic con- 
tinuous transition of rigid Ising-like systems (m=l) is fluctuation-driven to a 
first-order transition once the MEC is accounted for. On the other hand XY 
(m=2) and Heisenberg (m=3) systems are unaffected by the compressibility 
and the MEC is irrelevant [7-9] . This particular aspect is in accordance with 
the early results since the crossover exponent for the MEC is a, the rigid sys- 
tem specific heat exponent and the interaction is relevant for rigid systems 
with a divergent specific heat, i.e. a > 0. 
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Experimentally, the coupling is exhibited by the sensitivity of the crit- 
ical temperature to an applied external pressure. First-order and continuous 
transitions, in some cases separated by a tricritical point, were reported in 
different systems [12-16]. 

In the present article we shed a new light on the behavior of com- 
pressible magnets. We emphasize the drastic role played by quenched bond 
disorder, such as dilution (obtained by replacing some of the ions by non- 
magnetic ones). We find that for Ising-like systems the quenched disorder 
has a crucial effect in preventing the runaway to the first-order transition 
caused by the MEC in pure systems. This drastic effect of the disorder orig- 
inates from the so-called Harris criterion: Weak disorder is relevant for rigid 
systems with positive a (namely precisely the same systems that are affected 
by the MEC). These disordered rigid systems exhibit a new critical behavior 
associated with a random fixed-point for which a < 0. Hence one possibility 
is that a weak MEC will not affect their behavior (another possibility being 
that the MEC has the dominating effect: It drives the transition discontinu- 
ous which is then somewhat weakened by the disorder). In the next Section 
(2) the former heuristic conclusion is confirmed by a detailed RG analysis. In 
Sec. 3 the MET is applied and its results are compared with those obtained 
by the same method for pure compressible magnets. In the last Section (4) 
the important conclusions are summarized and discussed within the context 
of other bond-disorder effects in low-dimensional magnetic models. 
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2 Renormalization group analysis 

The momentum-shell RG analysis begins from a continuous-spin Landau- 

— * 

Ginzburg- Wilson Hamiltonian. With (f)(x) representing the m-component 
spin field and u{x) the d-dimensional local displacement the Hamiltonian is 

[4,7]: 

-(3H= J d^x{-r{xmx)f + -J2{VM^)) +Mm) + (1) 

a=l 

{^K - i/x)(V ■ u{x) f + /X X: (V«.(x))' + g{V ■ u){m f} ■ 

The first three terms represent the rigid random-bond Ising model. The 
dilution is accounted for in r[x) = r+6r{x), where r is the average of the 
local random "temperature" r{x). The local deviations obey 

< 5r{x)6r{y) >=< 5r^ > 6{x — y) , 

(higher moments are irrelevant and so are the omitted disorder-dependent 
terms in the elastic part of the Hamiltonian). K and /x are the bulk and shear 
moduli and g is the MEG. To average over the disorder we use the "n goes to 
zero" replica trick. After replicating and averaging the local displacements 
may be integrated out in < Z** >. The outcome is an effective replicated 
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Hamiltonian depending only the n m-component vector 0", a=l,2,....,n: 

r 1 ^ 1 ™ ^ 



where: 



U — U{i 



A{\K+{^)„} 



_ 1 2 1 2 



and 



A = 8 < > . 

The initial (and physically acceptable) values obey ti > 0, A > 0, v < 0. V 
is the volume of the system. 

The Hamiltonian has three quartic couplings. The flow within the 
three dimensional subspace of parameters which contains u, A, and v deter- 
mines the flxed-point structure of the system. Integrating out small scales 
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fluctuations, we obtain the recursion relations to order one loop: 



du 



m + 8 



■u^ + 2Au , 



(3) 
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mn + 8 



6 



m + 2 



3 



Au, 



(4) 




dv m 



m + 2 
3 



n + 2 
3 



(5) 



uv + 



One observation is called for at once: the coupling v cannot feed into 
the renormalization of the two others because of its special non-local form. 
Consequently we separate the search for the fixed points (f.p.) into two 
steps: We look first for the f.p. in the u-A directions. Secondly, with these 
f.p. values for u and A inserted into Eq. (5) we look for its fixed-points. The 
first step is analogous to the analysis done in the study of rigid random-bond 
models. The results are also m-dependent [17]. For the Ising model the pure 
f.p. is unstable towards the inclusion of the disorder coupling A, while it is 
stable for systems with m > 1 [17,18]. We therefore discuss the results for 
these cases separately. 
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2.1 Ising model (m=l) 

The stable f.p. in the u-A plane is the so-called Khmel'nitskii's f.p. [18] at 
which both u and A are of order ^fe. We therefore have to check the stability 
of this f.p. while switching on v. To that goal we insert the values of the 



two-loop calculations [18,19]. For v close to zero the term in is negligible 
and thus we have: 



Terms of order e were neglected compared with y^. It should be noticed 
that the crossover exponent for i> is a of the rigid random-bond system. This 
is contrary to the case in which a random-field is applied to the random-bond 
system and for which the crossover exponent is not 7 of the bond-disordered 
model [20]. Indeed, in the present case the v term is the square of the 
internal energy (per unit volume) and hence has a as its scaling exponent. 
Our finding dv/dl < implies that this fixed point is the stable one and 
hence is describing the critical behavior of this system (at least within the 
renormalized perturbation theory which starts from weak couplings). We 
thus conclude that elastic random-bond Ising systems will exhibit, at their 
critical point, the critical exponents of the rigid random-bond Ising system. 
If the MEG coupling v is large and that of the disorder A is small, the initial 




These values were obtained from 
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flow will be towards larger v. The stability boundary may be reached before 
the flow reverses itself to smaller values. In this case a first-order transition 
is still possible and the continuous transition will be attained (through a 
tricritical point) at a finite strength of the disorder. That is very different 
from the behavior in absence of disorder for which it was found that the 
elasticity is relevant at the order e Wilson-Fisher f.p. (of the pure-rigid model 
[7,8]) and drives the transition first-order. The flow diagram is schematically 
depicted in Fig. 1. 

2.2 Continuous spin models (m > 1): XY, Heisenberg, 
etc. 

For these models with continuous symmetry of the order parameter the stable 
f.p. in the u-A plane is that of the pure system. We therefore have to insert 
the values at this f.p. u* — [17] into Eq. (5): 



dv . 6(m-F2) , ,4-m, 

9Z ^ 3(m + 8) ^ V + 8^ ' ^ ' 

Although dv/dl is positive for m < 4 and small e at the u = f.p., it is again 
given by a of the rigid (and pure for the considered) systems. As discussed 
above this is a general relation which holds to all orders in e. It is known [17] 
that a becomes negative for larger values of e and is deflnitively so at e = 1 
for all m > 1. Therefore dv/dl is negative in 3d and the same f.p. of the 
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pure and rigid systems also describe their critical properties in presence of 
both disorder and compressibility (and their exponents remain unchanged). 



To go beyond the renormalized perturbation expansion we need to have re- 
course to mean-field theory (MFT), although its results arc always to some 
extent questionable (e.g. they are insensitive to the dimensionality of the 
embedding space). However, we may gain some insight by comparing within 
the same approach the predictions for the critical behavior with and without 
dilution. 

MFT is apphed to the lattice model with S'^ — 1 on every site. We 
choose a somewhat simplified (Bomb's) model which was shown [3, 5] to 
capture the essential physics. In this model all nearest-neighbor exchange 
interactions are strengthened (or weakened) by the same amount depending 
on the average inter-site distance. The deviation in this distance is propor- 
tional to the change in the volume fraction w — V/N [N is the number of 
sites) from its average wq. 

The Hamiltonian contains two contributions: the elastic and the mag- 
netic. The elastic part is: 



3 Mean Field Theory 




(8) 
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where 4> is proportional to the elastic constant. The magnetic part is [3]: 

Hmag = Ji,A^)SiSj , (9) 

where we have chosen Si to be a scalar (Ising spin) keeping in mind that MFT 
yields similar results for m-vector spins. The nearest-neighbor couplings are 
given by: 

Ji,j{w) = €ij{Jo + Ji{w - Wo)} , (10) 

ejj are independent random variables with the distribution p(e) = x6{e) + 
(1 — x)5{e — 1), where x is the concentration of non-diluted bonds (similar 
results will be obtained for either site or bond dilution). Ji is the coefficient 
of the part which depends on the inter-site distance (only the dominant, 
linear, dependence has been kept). 

It is straightforward to integrate out the volume fluctuations. The 
remaining effective Hamiltonian depends only on the spin degrees of freedom: 

- f3H,ff = 13 Jo E U,ASj + |^(E ^^,ASjf ■ (11) 

i,3 ^ i,j 

Introducing mean-fields, both terms are reduced to single-spins ones. Ex- 
panding the free energy and performing the average over the disorder yields 
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the following effective mean-field Hamiltonian: 



N J"^ N 

— (3Hmf = PJoxcm Si pjoxcrn^ + l3-\x'^c^m^ Si -/3 Jfx^c^m^ , 

Y 2 4(1) ^ 8(1) 

(12) 

where c is the coordination number of the lattice and here m is the magne- 
tization per site. ^Prom this Hamiltonian all thermodynamic properties can 
be calculated. In particular for the averaged free-energy we obtain: 



PfMF = lim = -^InZ (13) 

PJqXcw? pjfx'^c^m'^ , ,r^T pJix'^c^m^^ 

— ■ h In coshipjoxcm H } . 

2 80 40 



The critical temperature is given by that at which the coefficient vanishes: 
kTc = xcJq. The transition becomes discontinuous when the coefficient 
becomes negative. Therefore the tricritical point is at kTc = xcJq and Jf = 
-^gf-- We see first that within MFT the critical temperature depends only on 
the dilution with the same dependence as that of the rigid system. The value 
of Ji/Jo beyond which the transition turns first-order does depend on the 
dilution: The more diluted is the system the larger this parameter must be in 
order to change the order of the transition. Hence the regime of second-order 
is enlarged at the expense of the first-order one the more diluted the system 
is. The phase diagram is drawn in Fig. 2. This trend is consistent with the 
RG results for the Ising systems in which a first-order transition was turned 
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second order by the dilution. 

4 Conclusions 

To summarize, we have shown that disorder (e.g. dilution) has a crucial 
role in compressible Ising magnets: While pure systems are driven first- 
order by coupling spin fluctuations to the elastic degrees of freedom, none of 
that happens in disordered systems. Their critical properties are insensitive 
to the compressibility and continue to exhibit a random-bond continuous 
transition. This behavior is consistent with that of the specific-heat. Indeed, 
the Harris's [21] and other [17,22] arguments imply that a. of the pure system 
is also the crossover exponent for the dilution, which is therefore relevant for 
Ising systems (and irrelevant for m > 1 vector- models). 

Consequently, Chayes et al. [23] have shown that in any random system 
a is negative. Hence we reach the general conclusion that the compressibility 
cannot be relevant at a second-order transition of a random system. Our 
conclusions are also consistent with those of Refs. [24] . There it was shown 
that bond-disorder suppresses (in d = 3) any tricritical point (it eliminates it 
in d = 2 and therefore 2d compressible magnets undergo always a continuous 
transition in the presence of disorder). 

Although the crossover from one universality class to another is slow 
[a being a small exponent), compressibility and bond-disorder are present 
in many physical systems. Therefore it will be very worthwhile to observe 
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experimentally the effects of dilution in particular on the transition in com- 
pressible magnets with uniaxial anisotropy. 

Acknowledgments 

We are thankful to A. Aharony and D.J. Bergman for helpful comments. 
One of us (Y.S.) is grateful for the hospitality of the Groupe de Physique 
des Solides at Universite Paris 7. He also acknowledges the Donors of the 
Petroleum Research Fund administrated by the ACS, for partial support of 
this research. 



13 



REFERENCES 



[I] O.K. Rice, J. Chem. Phys. 22, 1535 (1954). 
[2] C. Domb, J. Chem. Phys. 25, 783 (1956). 

[3] D.C. Mattis and T.D. Schultz, Phys. Rev. 129, 175 (1963). 

[4] A.I. Larkin and S.A. Pikin, Zh. Eksp. Teor. Fiz. 56, 1664 (1969) [Sov. 

Phys. JETP 29, 891 (1969)]. 
[5] S. Sahnas, J. Phys. C7, C20 (1974). 

[6] V.B. Henriques and S. Sahnas, J. Phys. C20, 2415 (1987). 
[7] J. Sak, Phys. Rev. BIO, 3957 (1957); J. Bruno and J. Sak ibid 22, 3302 
(1980). 

[8] F.J. Wegner, J. Phys. C7, 2109 (1974). 

[9] D.J. Bergman and B.I. Halperin, Phys. Rev. B13, 2145 (1976). 
[10] A.F.S. Moreira and W. Figueiredo, Phys. Rev. B46, 2891 (1992). 

[II] See, e.g., Ref. [9] for a comprehensive overview. 

[12] C.W. Garland and R. Renard, J. Chem. Phys. 44, 1130 (1966); B.B. 

Weiner and C.W. Garland, ibid 56, 155 (1972); C.W. Garland and R.J. 

Pollina ibid 58, 5502 (1973). 
[13] C.W. Garland and B.B. Weiner, Phys. Rev. B3, 1634 (1971). 
[14] C. P. Slichter et al, Phys. Rev. B4, 907 (1971). 
[15] E. Erwin et al, Physica B161, 260 (1989). 

[16] H. Yoshida et al, J. Magn. Magn. Mater., 90 + 91, 177 (1990). 
[17] A. Aharony in Phase Transitions and Critical Phenomena, Vol. 6, 

14 



p. 358 (C. Domb and M.S. Green - Edit., Academic Press, N.Y. 1976). 
[18] D.E. Khmel'mtskii, Zh. Eksp. Teoret. Fiz. 68, 1960 (1975) [Sov. 

Phys. JETP 41, 981 (1976)]. 
[19] G. Grinstein and A. Luther, Phys. Rev. B13, 1329 (1976). 
[20] A. Aharony, Euro. Phys. Lett. 1, 617 (1986). 
[21] A.B. Harris, J. Phys. C7, 1671 (1974). 
[22] Y. Shapir and A. Aharony, J. Phys. C14, L905 (1981). 
[23] J.T. Chayes et al, Phys. Rev. Lett. 57, 2999 (1986). 
[24] M. Aizenmann and J. Wehr, Phys. Rev. Lett. 21, 2503 (1989); 

K. Hui and A.N. Berker, ibid, 2507 (1989). 



15 



FIGURES CAPTIONS 



Fig. 1. The three-dimensional flow diagram with the fixed points: 

(A) Pure-Rigid and (B) Disordered- Rigid. 
Fig. 2. Schematic mean- field phase diagram {Y = -^). 
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